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We discuss the description of full counting statistics in quantum transport with a nonMarkovian 
master equation. We focus on differences arising from whether charge is considered as a classical or 
a quantum degree of freedom. These differences manifest themselves in the inhomogeneous term of 
the master equation which describes initial correlations. We describe the influence on current and in 
particular, the finite-frequency shotnoise. We illustrate these ideas by studying transport through a 
quantum dot and give results that include both sequential and cotunneling processes. Importantly, 
the noise spectra derived from the classical description are essentially frequency-independent and 
all quantum noise effects are absent. These effects are fully recovered when charge is considered as 
a quantum degree of freedom. 

PACS numbers: 73.23.Hk, 73.23.-b, 73.63.Kv, 42.50.Lc 
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I. INTRODUCTION 

Stochastic processes are seldom Markovian: most 
stochastic models arise from reduction with respect to 
some large collection of degrees of freedom, the evolv- 
ing state of which does impart some memory of the sys- 
tem's past history. Despite this, most model systems 
are assumed to have Markovian dynamics, largely due to 
the inherent difficulty in treating nonMarkovian (NM) 
effects. If we restrict ourselves to the study of systems 
which can be described by a density matrix p(t) , ageneric 
nonMarkovian master equation (NMME) reads [l, [|J 



dt'W(t - t')p(t') + 7 (t). 



(1) 



This equation describes the evolution of p(t) under the 
action of NM kernel W(i) from a time to, and also in- 
cludes inhomogeneous term 7(i), which describes the ef- 
fects of the systems memory of its history prior to time 
to- For the important case in which Eq. ([1]) is obtained 
by tracing out bath degrees of freedom, 7(i) describes 
initial correlations between system and bath. These cor- 
relations typically decay in the long time limit and thus 
7(t —> oo) — >• 0. In a Markovian system (or Markovian 
approximation to the above), W(t — t') ~ S(t — t'), mem- 
ory effects arc neglected and the inhomogencity is absent. 

In this paper we focus on nonMarkovian effects in elec- 
tronic transport where Eq. (|TJ) might, for example, de- 
scribe the evolution of the density matrix of electrons 
in a quantum dot connected to leads. We are inter- 
ested in the statistics of the number of charges n trans- 
ferred through the system in time t — to, the so-called 
full counting statistics (FCS) [H, H[. This information 
is usually encapsulated in the moment generating func- 
tion (MGF), £/(X)i), whose derivatives give the moments 
(n«(t)> - i§#| x =o. 

Under a classical picture of electron transfer, Eq. (JTJ 



can be generalized to its n-resolved counterpart [H, Q 

A ft 

p (")(t) = V / dt'W {n - n \t-t')p^ n \t')+^ n \t), 
dt , ./*„ 



(2) 

which is an equation of motion for p( n \t), the density 
matrix of the system (e.g. QD) conditioned on n electron 
transfer events having occurred. The sum of all such 
partial density matrices is the full density matrix p(t) = 
J2n P (*)• I n tn i s equation, W^™~™ ' is that component 
of the kernel which transfers n — n! charges, and 7^(t) 
the n-resolved inhomogeneity, which describes the effects 
of initial conditions for each value of n. The probability 
of n charge transfers is obtained by tracing over system 
degrees of freedom: P(n; t) = Tr{p(' 1 ' (t)}, and the MGF 
is simply Q{x, t) = ^ n P(n;t)e mx . This approach is 
similar to the Markovian theory 0, H| . 

In reality, electron number is a quantum-mechanical 
variable and electrons can, for example, exist in super- 
position between states in the dot and states in the lead 
where counting takes place. The counting statistics of 
such a quantum degree of freedom was addressed by Lev- 
itov and coworkers by defining a MGF which yields the 
moments for a specific detection scheme Q , an approach 
which has since been widely applied and has been the 
subject of experimental interest pLQj | . 

The aim of this paper is to discuss how this quan- 
tum theory of FCS can be described with an n-rcsolvcd 
NMME of the form Eq. and to compare and con- 
trast quantum and classical analyses. In the quantum 
treatment, Eq. ([2]) arises as a reduced description from 
the tracing out of reservoir degrees of freedom from a 
fully quantum-mechanical system-reservoir theory. As 
we show, the difference between classical and quantum 
counting schemes is manifest in the n-resolved inhomo- 
geneity, j( n \t), which is uniquely determined in each 
case by calculating the MGF in the stationary limit. 
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This difference is attributed to the occurrence of system- 
reservoir superpositions at t = 0, obviously absent in 
the classical case. We also discuss the interpretation of 
p( n ^ (t) in the quantum case. 

We note in passing that the important role of initial 
system-bath correlations and memory effects is, of course, 
not limited to n-resolved NMMEs. An important exam- 
ple is the spin-boson model [Tlj , a prototypical model in 
the study of decoherence induced by external baths. 

Once the above points have been clarified, we discuss 
how the physics at finite frequencies is affected by the 
classical or quantum character of ^ n \t). In particular, 
we will focus on the current (I(t)) = 4f(n(t)} and its fluc- 
tuation spectrum S(w) = J™ x d,Te iut \{{8I(t), (5/(0)}), 
where SI(t) = I(t) — Id c measures deviations away from 
the steady-state current i^c = (I(t))\t-yoo, and we have 
symmetrized the definition [121 ] . This spectrum con- 
tains information about the internal dynamics of the 
conductor not accessible via dc measurements 
Importantly, we will demonstrate that the classical de- 
scription, only valid for large voltages where n-rcsolvcd 
Markovian QMEs arc justified @, H|, fails in describing 
frequency-resolved noise in situations where the measur- 
ing frequency is larger than both the temperature and 
the applied voltage, namely tuo > ksT, eV. In this 
regime, which has been already achieved experimentally 
\\H . quantum noise, i. e. zero-point fluctuations, domi- 
nates over thermal and shot noise. 

As example, we will study the noise spectrum of a 
quantum dot with a single Zeeman-split level in the 
strong Coulomb Blockade regime. In particular, we will 
consider the cotunneling regime using 4th-order Liouvil- 
lian perturbation theory (LPT) [L9l - |2ll ]. which allows us 
to go beyond the sequential noise results of e. g. Refs.JH- 
|2~H . We show that the noise calculated from quantum 
counting exhibits a series of steps not only at frequencies 
hoj = \ ± eV/2 — e a \, corresponding to transitions directly 
between the dot levels (energy e a ) and the leads, but also 
at ftw = inelastic with -Bi„ c iastic = \e± - Cfl tne inelas- 
tic spin-flip transition energy. Broadening of all these 
quantum noise processes by elastic cotunneling is also 
observed. In comparison the noise profiles derived from 
the classical description are extremely flat as functions of 
frequency and all quantum noise effects are absent. 

The paper is organized as follows: In Section|TT]we con- 
sider the classical NMME and derive an expression for 
the classical stationary inhomogeneity (Eq. (fl4)) ). This 
expression, which was obtained previously by Flindt et 
al. in Ref. docs not capture quantum noise effects. 
In Section IIIII where we explain how to obtain the quan- 
tum inhomogeneity (Eq. (|29|) ) by using Liouvillian per- 
turbation theory (technical details are relegated to the 
appendix). The current and noise spectrum obtained 
with our approach are described in Section IIVI and our 
results for the quantum dot example are described in 
Section |Vj Finally, Sec. IVII discusses the interpretation 
of the quantum NMME and the relation with previous 
results in the literature. 



II. CLASSICAL INHOMOGENEITY 

Classically the stationary MGF for the number of 
charges transferred to the collector in the interval to to 
tf is simply 

GM,t f o) = (e lx[nM - nM] )st,t (3) 

where n(t) is the collector charge at time t. At time to 
we require that the system be stationary. We want to 
calculate this quantity for a system whose dynamics are 
determined by a NMME. Let p( n > (t) be the partial system 
density matrix associated with presence of n electrons in 
the collector; the sum of all such is the full density matrix 
of the system: p(t) = p^ n \t). We specify the classical 
NM model by stating that if the system is prepared at 
some initial time ti in state p^ ni ^(ti) with collector 
charges, the subsequent evolution follows the n-resolvcd 
NMME 

h r* 

—p^(t) = J2 dt'W {n - n ] (t - t')p (n \t'). (4) 

There is no inhomogeneity here because we explicitly 
prepare the system without history. Introducing the 
counting field \ through the Fourier transform p(x, t) = 
E n e lnx p( n \t), we obtain the Laplace-transform solution 

/>oo 

p( X ,z)= / dt fl e- zt fp(x,t f )^n(x,z)p(x,t i ), (5) 
Jo 

with t xy = t x —t y and x _ dependent propagator f2(%, z) = 
[z — W(x, z)}^ 1 ■ At time tf > ti we have then 

P <n/) (t/)=E n ( n /*'*/> (n4) (*i)- (6) 

Tli 

with n xy = n x — n y and propagator Q(n, t) the inverse 
Fourier and Laplace transform of f2(x, z). The full sys- 
tem density matrix is p(z) = Y^ n P ( z ) = P(X = 
and, taking the long time limit, the stationary density 
matrix is p s tat = p(t — > oo) = linxj-^o zp(z) which sat- 
isfies W,o s tat = with W = W(x = 0;z = 0) and is 
assumed unique. 

To calculate the MGF of Eq. © (for arbitrary t f0 ) we 
require not the evolution from these prepared initial con- 
ditions, but rather from initial conditions corresponding 
to stationarity. This implies a history for the system and 
hence a finite inhomogeneity. To this end, let us consider 
the system evolving under Eq. Q at some intermediate 
time to with tf > to ^> ti. Were the system Markovian, 
we could use the semi-group property of the propagator 
f2(n, t) @, H[ to write Eq. (JB|) equivalcntly as 

P {nf) (t f )= n(n fo ,tf )n(noi,toi)p( ni Hu). 

no ,n,i 

With a NM kernel, however, this semi-group property 
disappears and we therefore write the evolution from ti 
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to tf via to as 

P (nf) (t f )= ^{nf ,tf )n(noi,toi)p {nt) (ti) 



ftf 

/ dt a n(n fa ,t fa )^(t toi, U 

„ Jtn 



(7) 



where the second term represents a "correction" to the 
Markovian result. In this case p^ n '{t) obeys the NMME 
of Eq. ^ starting at time t in state p^ n °^(t ) = 
Y, n . n(noi, toi)p (n 'Hti) with inhomogeneity 7 (n) (t, t oi , U) 
whose arguments make explicit its dependence on previ- 
ous times. 

To obtain an explicit expression for the inhomogeneity, 
we equate Eq. (j6]) and Eq. ([7]) and perform a Fourier 
transform of variable n/ to obtain 

tt(x,tfi)p(x,U) = tt(x,tfo)ft(x,toi)p(x,U) 



+ / <ft fi(x,*/a)7(x; *„, fa,**) (8) 

J to 

A Laplace transform over the two distinct time intervals, 
io°° Jo°° dtfodtoie-*^-^*, yields 



1 



{n( x ,z)-n(x,s)}p(x,t l ) 



z — z 

= n( X , z)Cl( X , S)p( X , U) + fi(x, z)j( X , z, z, U), (9) 
and rearranging, we obtain 

7( X , z, ~z, ti) = j—z (W(x, ^) - W(x, z)) a( x , z)p( X , u). 

Taking the long time limit of the to — U interval and iden- 
tifying the state of the system at to, the inhomogeneity 
reads 
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w=E r cr io W" o) (to) 



r ir ] ( z ) = - (w (n) (o) - c*)) . 

where the limit toi — > oo is understood here (and the 
state of the system at U becomes irrelevant). This result 
in Eq. ([7]) gives 



(t f )=J2^ f o,t f o)p^(t ) 

+ £ f 1 dt a ^{n fai t fa )T^ a \t a o)p^°\to). 



with the Fourier-transformed density matrix 

p(x;t/o) = E eixn/ V (n/) (t/) 



= ^ e ™|^Jl(n / o ) t/oy" (]) (to) 
+ z2 dt a n(n fa ,t fa ) 

n a ,no •' t ° 

x T^ \t a o)p {no) {to)\. (12) 



Note that the counting field x nere is associated just with 
the interval t fo and not t fi as in Eq. ((5J) . Finally, Laplace 
transform with respect to the interval t fo gives 



p{x,z) = £i(x,z) {i + r cl (x,z)}p s tat, 



(13) 



where in the last line we have used (t°) = 

p(to) = Pstat in the long-time limit. 

The n-resolved density matrix p( n > (t) that corresponds 
to Eq. (Jl3j) is the solution of an ra-resolved NMME like 
Eq. ^ with inhomogeneity 



7ci(x,z) = r c i(x;2 : )pstat; 
zr cl (x Z ) = W(x;0)-W(x;^). 



(14) 



Note that here n is the number of charges transferred in 
the interval between t and to (and not the total number 
of collector charges). The complete inhomogeneity is 



Tel 



W = E^r ) W = 7ci(x = 0;^) 



-W(z) 



Pstat(15) 



which disappears in the long time limit, lim4_ ) . 00 j(t) 
lim z ^ zj{z,t ) = -W(0)p st at = 0. 

Tracing Eq. (fl"3|) over n we obtain the MGF: 



Q C1 (X, z) = Tr {Sl( x , z) [1 + r cl ( X , z)} Pstat } . (16) 



With the inhomogeneity of Eq. (fT5|) and employing a no- 
tation in which ((...))= Tr {.. . p s tat} denotes the expec- 
tation value in the steady-state, the classical MGF reads 



Gc\.{x,z) = - + 1 



1 



z z z— W(x, z) 



W(X,0)». (17) 



The inhomogeneity of Eq. (fT4"]) was given by Flindt et 
al in @. and the classical MGF, like the inhomogeneity 
itself, can be expressed solely in terms of the x-resolved 
kernel. 



III. QUANTUM INHOMOGENEITY 



The MGF for the number of charges transferred in the 
interval tf to to is then 



The quantum- mechanical MGF of Levitov and cowork- 
ers M can be written as 1251 



Gd(x,tfo) = Tr{p(x,f/o)} , 



(11) 



^ q m(x,t)-(e- l *"e^ t ) e -^ ft ) s 



(18) 
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where n is the number operator of collector electrons. 
This form was originally derived by considering the pre- 
cession of a spin coupled to the current flowing through 
the device but any ideal passive charge measurement 
should give the same result. In this section we will con- 
sider this expression in the context of a system-reservoir 
theory and show how the MGF may be obtained from an 
n-resolved NMME of the form Eq. The inhomogene- 
ity in this case is very different to that found from the 
foregoing classical analysis. 

Let us consider a model Hamiltonian composed of 
reservoir, system, and interaction parts: H = H les + 
Hg + V. In its diagonal basis, the system part reads 
Hs = J2 a Ea\a)(a\, where \a) is a many-body system 
state of N a electrons. We assume noninteracting reser- 
voirs H rcs = Ysk^ka + fia)a ka aka, where a ka is anni- 
hilation operator for an electron of energy oj^a in lead a, 
and [i a is the chemical potential of lead a. Single-electron 
tunnelling between system and reservoirs is described by 
the Hamiltonian 

V = t kamala d m + 'L4 fl *oi ( 19 ) 
kam 

where d m is the annihilation operator for single-particle 
level m in the system, and tkam is a tunnelling amplitude. 

The density matrix g(t) of combined system and leads 
evolves according to the von Neumann equation: 

g(t) = -i[H,g(t)]=£g(t), (20) 

which defines the Liouvillian super-operator C. Eq. (|20p 
is solved simply as g(t) = IL(t — to)g(to) with full propa- 
gator II(t) = e ct . 

The MGF of Eq. jTSJ can be rewritten as 

Q{x,t) = {e iH ^ t e- iH ^% tat , (21) 
where the gauge-transformed Hamiltonian reads 

H( X ) = H S + H les + V( X ) (22) 
with counting fields appearing in the tunnel Hamiltonian 

V(x) = ^™ eiA,1 4A + H.c, (23) 

kam 

where we here assume that we count only in lead 7 (hence 
the Kronecker delta). The x- resorve d Liouvillian super- 
operator is defined via 

C(x)p = -i{H{\ X )p-pH{-\x)} (24) 

such that, with propagator II(x,i) = e £(x)t , the MGF of 
Eq. (JUJ) reads j24| 

<7qm(X>*/o) = (n(x,t/o)) s tat 

= lim Tr {n( X ; t/o)H(x = 0; t 0t )g(U)} (25) 

toi— >oo 

The MGF involves two time-evolutions: first without 
counting such that the complete system reaches steady 



state at time to, and then from time to to time tf 
with counting where we evolve the system with the x~ 
dependent propagator until time tf. For convenience 
we choose the initial state g(U) to be separable g(0) = 
po'X'Peq 3 w hh p the initial state of the system and p£ q s the 
initial state of the leads, assumed to be thermodynamic 
equilibrium. In Laplace space this MGF reads 

S qm (x;2)= lim z o Tr{n(x;z)n(0;z o )£(io)}. (26) 

We perform the trace over the lead degrees of freedom 
of Eq. (|26p using the technique of LPT — the propagators 
are expanded as power-series of Cv(x) an d contractions 
between the various tunnel vertices are then evaluated. 
Details of this calculation are given in appendix [XI but es- 
sentially there are three possibilities: contractions solely 
within the leftmost propagator give rise to a x-dependent 
reduced system propagator at frequency z, $!(%; z); con- 
tractions solely within the rightmost propagator give rise 
to the reduced system propagator at zero-frequency and 
without counting field; finally, contractions between the 
two original propagators give rise to a block which cannot 
be resummed and this is the inhomogeneous term. After 
tracing out the leads and taking the long-time limit, the 
MGF in Laplace space reads 

e qm (x; z) = m(x, z) (1 + r qm (x, «)))>, (27) 

which is clearly of the same form as the classical expres- 
sion Eq. (|17|) but with different individual terms. The 
system propagator is fi(Xj z ) = [ z — VV(Xj z )]~ l which in- 
cludes the x-dependent NM effective system Liouvillian 
W(x; z) = Cs + S(x; z) with self-energy S arising from 
the coupling to the leads. In terms of the LPT diagrams 
described in appendix |A1 to lowest order the self-energy 
reads (next-order diagrams are given in appendix) 

E( X ,z)=X^X + ..., (28) 

z 

where X denotes the system part of a counting (i. e. x~ 
dependent ) tunneling vertex; the horizontal line, a free 
propagator (with frequency argument 2); and the over- 
line denotes a contraction between reservoir operators. 
The stationary kernel without counting is W = W(x = 
0, z — 0) from which the stationary system density ma- 
trix is given by Wp sta t = 0. 

In the same language, the expansion (also to lowest 
order) of the inhomogeneity block reads 

zr qm ( X ; z) = (x^G) - + ■ ■ • ' ( 29 ) 

where G is the system part of a tunneling vertex without 
counting. It is thus clear that the diagrams of r qm (x; z) 
do indeed arise from bath contractions between two full 
propagations, one with and one without counting. This 
has to be contrasted with the self-energy of Eq. (|28|) 
in which all vertices are x-dependent. This emphasizes 
that the quantum inhomogeneity, and by extension the 
MGF cannot, in general, be couched solely in terms of 
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S(X, z). This is in contrast with the classical case where 
knowledge of the kernel W(x; z) is sufficient to determine 
the inhomogeneity and the MGF. 

The moment generating function can be written as the 
trace over the x-resolved density matrix 



p(x, z ) = fifo z ) t 1 + r qm(x, z )\ Ps 



(30) 



As in the classical case the corresponding n-resolvcd den- 
sity matrix is the solution of an NMME of the form 
Eq. ([2]) with inhomogeneity: 



1 qm stat ■ 



(31) 



We refer to the discussions section concerning the in- 
terpretation of these results in terms of an n-resolved 
NMME. Setting x = into expression Eq. (JSSJ gives the 
full quantum inhomogeneity to be 



lqm(z) = r qm (0; Z)p s 



-W(z)P*t 



(32) 



which is identical to the expression found in the classical 
case. 



IV. CURRENT AND NOISE 

We now derive expressions for the current and noise 
based on our two results for the MGF, classical (Eq. (jTT)) ) 
and quantum (Eq. (|27[1 ). The average stationary current 
is defined as the rate-of-change in the number of collector 
electrons: (I(t)) = ^j(n(i)} In Laplace space we have 



(/(*)> = z(n(z)) = z 







d(i X ) 



(33) 



The corresponding finite-frequency noise can be quanti- 
fied with correlation function 

r°° 1 

S(lj) = dte^-({6I(t),SI(0)}) (34) 

J — CO " 

with 51 = 1— (I) and where we have symmetrised the 
definition [l2j]. Using MacDonald's formula [28| this can 
be expressed as 



S(u) =lo dtsm(ujt)— [{n 2 (t)) - (n(t)} 2 ] (35) 

J — oo 

and, following Ref. as 



S(u) 



2 d(i X y 



{(Q(x, z = iu>)) + (z = -iu)} . 



In discussing results, we will consider the total noise, 
rather than that in a single contact. According to the 
Ramo-Shocklcy theorem, for a two-terminal conductor 
the total current can be written in terms of the currents 
through the left and right terminals as: /tot = o^Il+PIr- 
This change can be affected by including counting fields 
for both contacts (xl and X R f° r left and right leads 
respectively) and making the substitutions x t ot = Xl + 
Xaccum = Pxl - a\R- Differentiation of the 



MGF with respect to Xtot then give the total cumulants 
rather than those of a single lead Q . In the following we 
consider symmetric barriers and set a = ft = 1/2. 

We begin with the classical case and exapnd the kernel 

as 



with 



W(x,z) 



W (k \z) 



k=0 



Qk 

d(ix) 



M 



:W( X ,z) 



(36) 



(37) 



The component W^°\z) = W(z) is the non-counting ker- 
nel and Q(z) = [z — W^ '(z)] is the standard propa- 
gator without counting field. 

Differentiating the classical MGF of Eq. (|17|) . we ob- 
tain the current (I} c iX z ) = z((W (1) (0))), such that the 
current 



(I) cl .W = (/)ci. = ((W«(0))), 



(38) 



is stationary and depends only on the first derivative 
block at zero-frequency Similarly, the noise reads 

Scl(w) = ((w (2) (o))} + «wWMttMwW(o))) 

+((w (1) HwM-^)w (1) (o)». (39) 

The first term is the autocorrelation term and is eval- 
uated at zero frequency. The remaining terms contain 
first-derivative blocks evaluated at both zero and finite- 
frequency. 

Taking the first derivative of the quantum MGF of 
Eq. (|27p . rearranging and canceling, we find the station- 
ary current 

(I) qm . = ((J (1 H;0))). (40) 
Here the current block (to sequential order) is 

jW(zi,zo) = G-X' + X' 



i 

G 



*0 



(41) 



where X' is the derivative of the system part of the X - 
dependent vertex evaluated at x — 0- As it appears in 
Eq. the leftmost frequency argument of this block 

is replaced with a dot since all diagrams with leftmost 
G-superoperator, and thus all those containing zi, are 
exactly zero in the stationary expectation value [29(. 

Similarly, taking the second derivative and simplifying 
we find the quantum expression for the noise reads 



+ (w -> -u>), 



(42) 



with J"' 1 ) block as as in Eq. (|4"TT) and with second-order 
block 



j( 2 )(z 2 ,Z!,Zo) 



2X 1 - X 1 

21 



X" - G + G - X" 

20 22 
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Expression Eg. (1421) reproduces the sequential expres- 
sions of Refs. U24J and, as is made explicit in the 



£0/k„T 



to/k„T 



appendix, extends them to cotunneling and in principle 
beyond. 

It is clear that the classical and quantum expressions 
for both the current and the noise have the same general 
block structure. The blocks themselves differ between 
the two approaches, however. In the classical case, one 
has the simple W-derivative blocks, whereas in quantum 
case, the J'-blocks are not directly related to the orig- 
inal kernel and must be calculated separately. We can 
directly compare classical and quantum results if we con- 
sider the classical Liouvillian W(x; z) to be that arising 
from our microscopic quantum model — the only differ- 
ence between the two calculations is then the nature of 
the counting and hence the inhomogencity. In this case 
the blocks W^ k ' and J 1 ^ have the same structure, but 
differ in the assignment of free-propagator frequencies 
within the diagrams. For example, the first-derivative 
block W^(z) has the diagrammatic expansion 



W {1) {z) = G-X' + X' - G + 



(43) 



This is the same as the block, Eq. (|4"Tj) but for the 
frequency arguments. 

Since, in both expressions for the stationary current, 
zero is the only relevant frequency argument, the two 
currents are equal 



Wcl = (i)c 



(44) 



This is as we expect since memory effects should be irrel- 
evant for the stationary properties of the system. In fact, 
since the inhomogeneity vanishes in the z —> limit, the 
zero-frequency FCS of classical and quantum approaches 
must be identical. We find the zero-frequency limit of 
the noise 



Sqm.(0) 



Sk(0) 

((W (2) »+2((W (1) RW (1) » 
+2(1) f((W (1) )) + ((W (1) RVV)) 



(45) 



with all blocks evaluated as zero-frequency, X = dX/dz. 
and R(z) the projection of the propagator Vt{z) out of 
the stationary state, which (up to definitions) agrees 
Ref. [H, HO]. At finite frequency, however, the MGFs and 
the noise results from the two approaches do differ. We 
explore these differences by considering an example from 
quantum transport in the next section. 



V. QUANTUM TRANSPORT RESULTS 

We consider the transport through a quantum dot with 
a single Zeeman-split level that we model with the An- 
derson Hamiltonian 



H 



7 d\d a + U n^n x + H rcs + V, 



(46) 
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FIG. 1. The finite- frequency total noise S(uj) as a func- 
tion of w for the infinite-?/ Anderson model QD with level 
positions e-f = — 15fcsT, ej_ = 5fcsT and chemical poten- 
tials fiL = -/in = \eV\ V/k B T = 2,15,25,50. Results 
are shown for the noise calculated with both quantum and 
classical counting from fourth-order (cotunneling) LPT ker- 
nels (solid black and dashed blue lines respectively). Both 
counting theories give the same noise at u = but as fre- 
quency increases, differences start to develop. This is most 
evident at small biases where, whereas the classical noise is 
essentially flat as a function of frequency, the quantum result 
shows marked steps at oj = \/j, a — ef\ (steps at ui — |/x a — 54, | 
are also visible at V = 50fcsT; see text). For V = 2ksT, the 
difference between between the two results is approximately 
two orders-of-magnitude at high frequencies. At higher bias 
the system enters the shotnoise regime and the relative impor- 
tance of the quantum contributions lessens. Other parameters 
were Tl = Tr = \k B T and bandwidth X c = 300fc s T. 



with e CT the energy of a spin- a electron in the dot and 
U the interaction energy. The reservoir and interac- 
tion Hamiltonians are as section III, with index a in- 
cluding both lead (= L 7 R) and spin (=t,4) index. For 
simplicity we take the limit of large Coulomb repulsion, 
U — > 00, and thus restrict the Hilbcrt space of the QD to 
three states: "empty" , and spin-up and spin-down single- 
electron states. We consider a bias symmetric about zero, 
Hl = —fJ-R = \cV , and two different level configurations. 
In both cases, we discuss the total noise S(u). 

In the first configuration we set tf = — 15ksT, = 
5ksT and sweep bias and frequency. Results for this con- 
figuration are shown in Fig. Q] and Fig. [2] as calculated 
using 4th-order LPT kernels (which include sequential 
and cotunneling contributions) and, in the case of Fig. [21 
quantum results calculated with 2nd-order LPT kernels 
(sequential processes only). At zero frequency, classical 
and quantum results coincide as expected. At finite fre- 
quencies, the contrast is stark. Concentrating on Fig. [1] 
first, we observe that, for low bias, the classical expres- 
sion predicts a noise that is essentially flat as a function 
of frequency. In contrast, the quantum noise increases 
in a step-wise fashion. These steps are the quantum- 
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FIG. 2. As Fig. [T] but with finite-frequency noise S(u) plot- 
ted as a function of applied bias for u/ksT = 0,4,10,23. 
Quantum- and classical-counting results with fourth-order 
LPT (cotunneling) kernels (black solid and blue dashed, re- 
spectively) are show alongside results from quantum count- 
ing with second-order LPT (sequential) kernels (green dotted 
lines). The quantum results are dominated by sequential pro- 
cesses with shotnoise step at V = 30fcsT and quantum noise 
step up at V — 30/cbT — 2u (u> < 15fcsT) and step down 
at V = — 30fcsT + 2uj for high frequencies (u) > 15k bT) in- 
stead. With the parameters cotunneling predominantly serves 
to increase noise, both at zero and finite frequency, at low fre- 
quencies. In contrast to the quantum case, the classical results 
change very little with increasing frequency 



noise steps described in Ref. [24j arising from quantum 
fluctuations between system and reservoir. With sequen- 
tial processes dominant, we expect quantum noise steps 
at uj = \fi a — e a \ for some a = L,R and a =t,l- For 

V < 30/cbT only the transition from the e-j- is visible 
since, due to the Coulomb blockade, this is the only level 
populated in the steady state. For V = 50fcsT both 
levels lie within the transport window and all four tran- 
sitions are visible. 

Fig.0sh.ows the noise as a function of bias. At a bias of 

V = 2e-|- the lower-lying level enters the bias window re- 
sulting in a step as we enter the shotnoise regime. Noise 
results at and above this step are insensitive to calcu- 
lational details since sequential electron hopping, which 
is essentially classical in the weakly coupled limit, then 
dominates. Below this voltage, quantum noise effects are 
visible. The dominant quantum noise process here is the 
excitation of an electron from the lower state (energy e-j-) 
to the leads. This is a sequential process and can occur 
for a bias V < 2(e T + uj) to the left lead (fi L = +V/2), 
and for V > — 2(e T + u>) to the right lead (/xr = — V/2). 
With V > and for w < |e T |, we therefore obtain a sin- 
gle step up at V = — 2(e^ +uj). For uj > |e T |, exitation 
to both leads is possible at V = 0, and low bias noise 
is correspondingly high. At V = 2(e^ + uj) one of these 
channels closes and we obtain a step down. 

Considering now both sequential and cotunneling pro- 



cesses together, the zero-frequency curve of Fig. [2] repro- 
duces the results of Ref. [2(| and, of course, classical- 
and quantum-counting results agree. The most signifi- 
cant new feature that occurs at cotunneling order is a 
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we asso- 



broad feature at low bias — following Ref. 
ciate this with inelastic cotunneling centered about the 
spin- flip transition energy -Emciastic = | e 4 — e f|- For these 
parameters, the main extra contribution of cotunneling 
processes at finite frequencies is to increase the fluctua- 
tions at very small bias. Quantum noise features directly 
attributable to quantum cotunneling are masked by se- 
quential processes here. The classical results at finite 
frequency differ strongly from the quantum case — at 
finite frequencies, the basic form of the classical noise is 
broadly the same as at zero frequency (sequential step at 

V = 2\ef\, flat inelastic cotunneling feature at low bias). 
There is a slight shift downwards of the classical curve 
for uj > compared with uj = 0, but this is a relatively 
small effect. 

The finite- frequency Fano factor F(uj) = S(u>)/(I) fur- 
ther serves to exaggerate the effects of quantum noise 
(not shown). At zero frequency, the two main features 
of the Fano factor as a function of voltage are a diver- 
gence of the Fano factor for eV — > 0, consistent with the 
fluctuation-dissipation theorem, and a pronounced super- 
Poissonian region strongly influenced by inelastic spin- 
flip cotunneling processes. 0, |3l|. In the classical case 
finite-frequency measurement serves to reduce the Fano 
factor such that noise becomes subPoissonian everywhere 
except for the divergence at V — > 0. In the quantum case, 
however, the superPoisonian peak increases in height and 
moves to lower biases. The result is massively increased 
Fano factor (F(u>) > 100) for a significant bias range 

V <j -^inelastic- 

Figure [3] shows results for a configuration in which 
both dot levels are far below the transport window (e T = 
— 20fcs7\ = — SObsT) such that d.c. transport can 
only proceed by cotunneling processes. The differential 
conductance G = d(I) /d(V) is plotted in the inset which 
would be vanishingly small with just sequential terms in 
the kernel, but is finite with cotunneling terms included 
here. Of particular note is a step in the conductance at 
the inelastic cotunneling energy V ~ inelastic [13 • 

The main panel of Fig. [3] shows the derivative of the 
quantum shotnoise dS qm (uj) / duj as a function of uj. Re- 
sults arc plotted from both second- and fourth-order ker- 
nels, labelled LPT2 and LPT4 respectively. The dom- 
inant feature in this plot is a large peak centered at a 
frequency u resonant with the transition between the 
lowest dot state and the chemical potentials, i.e. for 
uj = \n a — ej. Since we are here near linear response, 
A*l ~ Mi?, we have a single peak — increasing the bias 
splits this peak in two (not shown). 

As comparison of the LPT2 and LPT4 results show, 
cotunneling processes modify this sequential peak in two 
ways. Elastic cotunneling processes serve to broaden 
the sequential resonance. More importantly, cotunneling 
leads to a significant increase in dS(uj)/duj on the low- 
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FIG. 3. Transport through the infinite- U Anderson model 
with both dots levels well below transport window: e-j- — 
— 20fc_gT, — — SOksT; other parameters as in Fig. [1] The 
inset shows the differential conductance G = d{I) /d(V) as a 
function of bias. A clear jump is observed around the inelastic 
cotunneling energy .Einelastic — \ € i — e tl — lOfcsT. The main 
panel compares the derivative of the shotnoise, dS{u)/duj, as 
a function of w for both sequential (LPT2) and cotunneling 
(LPT4) kernels near linear response (V = lfcsT). Both pre- 
dict a large peak when the measurement frequency is approxi- 
mately equal to f^L/R — ^i ~ 30A;bT. In the LPT4 results, this 
peak is broadened by elastic cotunneling processes, but more 
importantly, there is a significant enhancement in the low- 
bias tail which can be attributed to resonance with inelastic 
cotunneling process u) ~ -©inelastic 



frequency side of the main resonance. This we attribute 
to the effects of inelastic cotunneling and this feature is 
situated around resonance with the inelastic cotunnel- 
ing energy u> = -Emciastic- We note that this feature is 
rather broad and does not have a simple peak structure. 
For comparison, the classical dS(ui)/dui is several orders 
of magnitude smaller than the quantum result (too small 
to plot on the same scale as Fig. [3]), and is approximately 
flat as a function of uj and shows none of these resonant 
features. 



ular the zero-frequency FCS, since 7(x, t — > oo) = 0. 
However, as we have discussed in general and for our 
quantum dot example, this difference is visible in the 
finite-frequency FCS and the influence on high-frequency 
quantum noise can be dramatic. 

The interpretation of the n-resolved NMME and the 
associated probability P(n,t) in the quantum case is 
interesting. As is clear from Eq. (|29|) . the inhomo- 
geneity contains counting field factors e tx ^ 2k+1 ^ 2 ^ ; k = 
0, ±1,±2, ... due to the pairing of counting (X) and 
non-counting (G) vertices (in the self-energy, counting 
vertices are always paired with other counting vertices). 
The presence of non-integer powers of e zx means that for 
the n-resolved NMME to contain the same information 
as its v-resolved counterpart, the transformation from x 
to n representation should be a continuous Fourier trans- 
form (cf. the classical case, where the n is clearly integer 
and the discrete transform is appropriate) and the above 
counting field factors imply an interpretation in terms of 
half-integer charge. Furthermore, in the quantum case 
there is no guarantee that P(n,t) is positive; The best 
interpretation of P(n,t) and, by extension, the diagonal 
elements of p^(t), is as quasi-probabiliti es l35l . |36| . 

Interestingly, Shclankov and Rammer [33f introduced 
the MGF 

G S R(x,t) = £<4e iH( -* /2)t e-^*/ 2 >'P n ), (47) 

71 

with P n being a projection operator onto states of defi- 
nite charge n (a similar MGF was considered by Esposito 
et al. [34( | ) that removes the interpretational problems 
of the MGF in Eq. (f2"Tj) . This is achieved by removing 
all traces of charge superposition in the initial density 
matrix. At sequential order, this projection approach 
is equivalent to our classical approach; at higher orders 
the approaches differ since Eq. (|4"T)) still admits super- 
positions within a charge-sector, which are absent in our 
classical approach. Nevertheless, since such superposi- 
tions do not contribute to the charge fluctuations at to, 
we anticipate very similar results for Eq. (|47p and our 
classical approach. 
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VI. DISCUSSION 

We have discussed how the description of the FCS of a 
transport system in terms of an n-resolved NMME differs 
depending on whether the number of transferred elec- 
trons is considered a classical or a quantum variable. 
The primary difference manifest is the inhomogeneity 
- in the quantum case, the inhomogeneity can not be 
described in terms of the system Liouvillian alone but 
rather requires new diagrams which reflect initial super- 
positions between system and reservoir that are obviously 
absent in the classical case. This difference is not observ- 
able in the long-time properties of the system, in partic- 
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Appendix A: Liouvillian Perturbation theory 

For completeness, we review here the basics of LPT 
and the translation of the resulting kernels into diagrams. 



More details are to be found in Rcfs. 
the diagrams. 
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To ease book-keeping, we introduce a compact single 
index "1" to denote the triple of bath indices (£1, k\, a\). 
The first index £1 = ± indicates whether a reservoir op- 
erator is a creation or annihilation operator: 

01 = a ilkiai = { f = + - (Al) 

L ak iai , ?i = — 

Leaving sums implicit, the reservoir Hamiltonian reads 

H rcs = (Uka + Pa) tt+tod-to 

= (wi + hi) aia T %,+ , (A2) 

where 1 denotes (— £i, k\, a±). In equilibrium, the reser- 
voir electrons are distributed according to the Fermi func- 
tion f(uj) = (e^/ k BT _|_ ]\ ^ which, since we include the 
chemical potential in Eq. (|A2[) and assume a uniform 
temperature, is the same for all reservoirs. 

We then define the system operators g^ a = 
y] m tkamjm ; such that the interaction Hamiltonian of 
Eq. (|19p can be written V = ^\a\gi. Correspondingly, 
the tunnel Liouvillian reads 



(A3) 



[V,-] = -itiiY J pcT p AlGl 



with bath-supcroperator A defined as 

and system supcropcrator G via 

G?0 = ^x( p=+ . (A5) 

The object a p is a dot-space supcropcrator with matrix 
elements 

, „> r r f 1, N s — N s > = even . . „. 

(o»).„ iir = { /( . ^ _ £ = odd , (A6) 

where, A, is the number of electrons in state s. 

The von Neumann equation reads g = Cg, with Liou- 
villian super-operator 

C = C Ies + C S + C V (A7) 

with £ rcs = -i[H res ,»], C s = -i[H s ,»], and C v = 
— i[V, • ]. The solution in Laplace space is g(z) = 
[z — -Cp 1 g{0). The reduced density matrix of the dot 
is given by tracing out the electron reservoirs 

Ps(z) = Tr R { p{z)) = Tr R | ^P(O)} • (A8) 

This we expand in powers of Cv to obtain 

p s (z) = Tr R {(n (z) + n (z)£yn (z) + . . .) p(0)JA9) 

with free propagator Hq(z) = [z — C rcs — Cs] 1 • With 
substitution of Eq. (|A3[) , a typical term of the expansion 
Eq. (|A"9j) reads 

H)" (ft^) TrR { noi^K^n ■ ■ • 

...a* A? Gf n o (z) P (0) }.(A10) 



Evaluating the action of the cr p superoperators we obtain 
a factor n° dd -P' an d, as the G operators also contain c, 
they evaluate at different positions in the chain as 



Gf'0=M l+1 x\ 9 J ) °> Pl = 
1 I Og u Pi = 



(All) 



We then bring all the A operators to the right of the G 
operators by using the relations A\G\, = —pp'G\,A\ , 
Tr R £ ros = , A-csPros = and A{C rcs = (C Tes - Xi)A{ 
with 

xi = -i£,i(uJi + /Lt Ql ). (A12) 
and our typical term becomes 

(odd \ 
Y[p, J Tr R { fl (z)G%> fV^-OG^Y . . . 

. . . G p 2 2 Qo{zi)G p 1 1 n (z)A^ . . . Af p(0)} (A13) 
with free dot propagator 

1 



fio(z) 



z - C s 



(A14) 



and z rn = z + J2i= m +i x ll 1 < m < n - 1. 

The the reservoir expectation values, 
Tr res {^...Af/^ s } - (^"...Af) eq) can be 
evaluated with the rules of Wick's theorem in Liouvillc 
space, which read [l9| 

• decompose (. . .) ros into pair-contractions 

• add minus sign for each interchange of A 

• omit factor ^Il°' dd P'') arising from a super- 
operators 

• each pair contraction then contributes a factor 

{AfAf) = 7 ?f 1 = <W>i/(-6Pi^i). (A15) 

With these rules, our typical term becomes 
(- l )"fi (^)G^^o(^-i)G5:'r 1 1 ■ ■ ■ 

...Gfn Q { Zl )Gfn Q {z) PS {t Q ) 

*( E (-1)^11 ^).(A16) 

\decomps / 

where the last factor indicates a sum over all pair de- 
compositions with the relevant Wick sign ( — 1)^. This 
expansion for p(z) can then be resummed to yield the 
familiar form 



1 



z- Cs - 



p(to) 



(A17) 
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with the self-energy 

even 



n \irrcd. / 

xG^n (z n . 1 )G p n n S 1 1 ■ ..G^n^ci 1 , 

where the sum is over irreducible contractions only. Typ- 
ically one approximates the full kernel by truncating the 
self energy at a given n. Here we consider only second 
and fourth order contributions. 

At second order, there is only one contraction which 
reads 



S( 2 )(^) = (-l)GfOo(^)G? 1 72T 1 - 
This translates into diagrams as 



S (2) (z) = G-G. 



(A18) 



(A19) 



which depicts two system-parts of the tunnel-vertices 
linked with a bath contraction. All forefactors are all 
left implicit in these diagrams. 

At fourth order, there are two irreducible contraction 
schemes: (41)(32) and (42)(31), which we label "D" for 
direct and "X" for exchange. The direct contribution 
reads 



£( 4D >(z) = Gfn (z 3 )Gfn (z 2 ) 

xGrO (z 1 )Gf 74 p f 1 7ff 2 - 
and the exchange term is 

£(«)( 2 ) = -Gi i n { Z3 )Gfn (z 2 ) 

Diagrammatically, we have 



(A20) 



(A21) 



r 



Y< {AD) {z) = G-G-G-G 
^ 4X \z) 



i 1 1 I 

G-G-G-G 



(A22) 



The evaluation of the various integrals required to imple- 
ment these kernels is discussed in Ref. [201 ] - 

Counting is effected by adding counting fields into the 
tunnel Hamiltonian such that, in our compact notation, 
we have V(px) = CiOiSie 1 ^ 1 */ 2 . Correspondingly, the 
X-dependcnt tunnel Liouvilian reads 

p 

= -i^p* p AlX?( X ) (A23) 
p 

which defines X P (x) = G\e lpSa ^ lXa ^ I 2 as the symbol for 
the system-part of of the x-dependent vertex. 



Performing the same expansion as before but with the 
X-dependent Liouvillian gives us the x~dependent non- 
Markovian self-energy to fourth order 



i 1 i 1 

X-X+X-X-X-X 



+X - X - X - X 



(A24) 



Appendix B: Quantum Inhomogeneity 

Writing the quantum MGF of Eq. (gBJ as Q(x, z ) 
Tr z)} and expanding the propagators we have 

p(x, z ) = lim zoTr R {n(x; z)II(0; z o )g(0)} 

zq^O 

= lim z Tr R {(IIo(z) + U Q (z)C x U (z) + . . .) 



x (U (z ) + n (z )C v U (zo) 



We proceed as in the previous section by evaluating con- 
tractions between the various tunnel vertices. There are 
three possibilities in doing this. Contractions can be 
made which gives diagrams which originate solely within 
the leftmost propagator. The vertices in such diagrams 
are all counting vertices and taking all such diagrams 
together, we obtain the reduced system ^-resolved prop- 
agator at finite frequency z. Similarly, diagrams with 
contractions solely from the rightmost propagator give 
rise to the plain (i.e. non-counting) system propagator 
at frequency zq — > 0. Finally, contractions can also be 
made between the two original propagators. The sum of 
all such diagrams can not be resummed and is thus left 
as a single block and this is what gives rise to the inho- 
mogencous term. Taking all diagrams into account, the 
final form of the MGF is 



g( X ;z) = hmz o ((n(x,z){l + T vn ( X ,z))n(0,z o ))) 



({n(x,z)(l + T qm ( X ,z) 



(Bl) 



Explicitly writing out only sequential and dircct- 
cotunneling terms, the form for r qm that one arrives at 
is 



1 



r qm (x; z) = x--g + x-x-x--g 

z z z z 



+X-X--G-G+X--G-G-G+ 

zzOO zOOO 

The exchange-cotunneling diagrams have the same fre- 
quency assignments as the three direct diagrams, but 
with the different pattern of contraction. The inconve- 
nient thing about this expression is the double free prop- 
agators. These can be removed by going one step back- 
wards and noticing that the product of free systcm-bath 
propagators can be split as 



n (z)n (z ) 



i 



z - ZQ 



(IIo(zo) - n (z)) . (B2) 
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Performing this separation and re-doing the contractions, 
we see that we can make the diagrammatic replacement 



2 




The inhomogeneity block becomes 



i 1 i 1 i 1 i 1 

zT qm ( X ; z)=X -G+ {-1)X -G + X -X -X -G+ (-1)X -X-X-G 



z z z 



I F=i 1 I F=i 1 I 1 I FF 1 

+X-X-G-G + {-1)X -X-G-G + X-G-G-G + (-l)X -G-G-G 

z z z 200 

Comparison of this form with Eq. (|A24[) allows us to write 



zr qm ( X ; z) = -E(x; z) + X-G + X-A 

2 



1 I FF I I F=. I r 



+X-G-G-G + X- A- G- G + X- X- A- G + X- X- X- A + ... (B4) 

000 200 220 222 

where we have introduced the shorthand 

A = X + (-l)G. (B5) 

This final form for the inhomogeneity is rather useful in evaluating the cumulants, because any A that remains 
undifferentiated gives zero when the x ~^ limit is taken. 

Appendix C: Current blocks 

Writing out explicitly only the sequential and direct-cotunneling contributions, the first and second current blocks 
read 



J (1 \ Zl ,z Q ) = G-X' + X'-G 



I r=^ I I i=i 1 I F=i 1 I FF 1 

-X'-G-G-G + G-X'-G-G + G-G-X'-G + G-G-G-X' + ... (CI) 

ZQ Zq Zq Z\ Zq Zq Z\ Z\ Zq Z\ Z\ Z\ 



1 (2) r t '/ 1 ^ ' 1 ' '// 

2^271 ( z 2i z l i z o) = X — X + 2~X" — G + ^ G — X 

Z\ Zq Z2 



I i=i I I F=. I I r=> I 

-X' -X'-G-G + X'-G-X-G + X'-G-G-X 

Z\ Zq Zq Z\ Z\ Zq Z\ Z\ Z\ 

I i=i I I i=i I I FF I 

-G — X' — X' — G + G — X' — G — X' + G — G - X' - X' 

Z2 Zi Zq Z2 Z\ Z\ Z2 Z2 Z\ 

I i=i I I 1 I FF 1 I FF 1 

-\X" -G-G-G + \G-X"-G-G + \G-G-X"-G + \G-G-G-X" + . .(.C2) 

2Q 2q 2o 22 2o 2o 22 22 Zq 22 22 22 
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